Fermions On One Or Fewer Kinks 
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We find the full spectrum of fermion bound states on a Z2 kink. In addition to the zero mode, 
there are int[2m//m s ] bound states, where m/ is the fermion and m s the scalar mass. We also 
study fermion modes on the background of a well-separated kink-antikink pair. Using a variational 
argument, we prove that there is at least one bound state in this background, and that the energy 
of this bound state goes to zero with increasing kink-antikink separation, 2L, and faster than e~ a2L 
where a = min(m s ,2m/). By numerical evaluation, we find some of the low lying bound states 
explicitly. 



I. INTRODUCTION 

A novel feature of fermion-topological defect interac- 
tions is the appearance of fermion zero modes [l|, [U, Q • 
The existence of zero modes has important implications, 
leading to phenomena such as fractional quantum num- 
bers [ij and superconducting cosmic strings ||. In any 
physical setting, however, the system is expected to con- 
tain both defects and antidcfects, and extended topolog- 
ical defects will frequently occur as closed structures, for 
example, closed loops of cosmic string, or closed branes 
in brane cosmology. Then it is important to determine 
the fate of a fermion zero mode in these situations. 

The fate of fermion zero modes on topologically trivial 
structures, such as kink-antikink or cosmic string loop, 
has been addressed in Ref. [f|. The expectation that 
the fermion zero modes would be recovered as the kink- 
antikink separation, or the size of the cosmic string loop, 
is increased indefinitely, was not met in Ref. Q. In the 
present paper, our primary aim is to reconsider the prob- 
lem of fermions on kink-antikink backgrounds. Contrary 
to Ref. we find that there are bound states on kink- 
antikink pairs whose energy vanishes exponentially fast 
with separation of the kink and antikink. 

We start by finding all fermion bound states on a sin- 
gle kink. If 2m f < m s where m / and m s are the fermion 
and scalar masses, we find that the bound state spec- 
trum only contains a zero mode. However, as we in- 
crease the fermion mass further, the number of bound 
states increases and is bounded by 2m f / m s as described 
in Sec. IIIIl We then turn to the kink-antikink system, 
proving first that a bound state exists if the kink and 
antikink are well-separated. Our proof is based on a 
variational argument and allows us to obtain an upper 
bound on the energy of the bound state. The bound it- 
self shows that the energy goes to zero with separation 
(2L) faster than exp(— a2L) where a = min(m s ,2m/). 
Next, we evaluate the bound state energies numerically 
and confirm the exponential dependence on L. We also 
find an exponential decay of the ground state energy with 
increasing 2mf/m s . 

In the next section we set up the problem. We sum- 
marize our results in Sec. [V] Identities involving hyper- 



geometric function are included in the Appendix. 



II. SETUP 

The 1+1 dimensional field theory we are interested in 
is described by the Lagrangian 
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- V 'Y + l^d aV J - gt/rW (1) 



where is a real scalar field, vj is a two-component spinor, 
and the 7' 1 are defined as 
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There are two masses in the model. The scalar mass is 
m s = \/2\ r q and the fermion mass is m/ = grj, where we 
are taking g > 0. 

The Zi kink solution has the well-known form (e.g. see 
Ref. 3) 



r/tanh (^y^) 



(3) 



and the antikink is obtained simply by letting z — > — z. 
We shall also be interested in the system that contains 
a well-separated kink and antikink, for which the scalar 
field configuration can be chosen to be 



77 tanh 



(-f(z + L) 



i] tanh [~^~(. z 



L 



n (4) 



The kink-antikink separation is 2L. 

Fermionic modes are found in the fixed scalar field 
background by solving the Dirac equation, 



(5) 



where we will consider (j> to be the kink solution of Eq. (|3|) 
and the kink-antikink configuration in Eq. The 
modes will contain a set of bound states (\E\ < m/) 
and continuum states. In this paper, we will only be 
interested in determining the bound states with E > 0. 
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Before proceeding further, it is convenient to perform 
a change to dimensionless variables defined by 
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2m t 



E' 



2E 



In what follows, we will drop the primes for notational 
convenience. The Dirac equations are then still given by 
Eqs. ([7]), (|H|), though with all variables having their di- 
mensionless meanings, and the (rescaled) kink and kink- 
antikink backgrounds read 

0k = tanh z (9) 



7 KK 



ee tanh(z + L) - tanh(z - L) - 1 (10) 



By substitution of one of Eqs. ([7]), ([8]) into the other, 
we obtain the 1-dimensional Schrodinger equations for 

-d 2 z f3±+g(g^Td z <P)f3 ± =E 2 p ± , (11) 

allowing us to identify the potentials 

V±(4>) = g(g4> 2 t d z <t>) (12) 

Note that Eq. (TlTj) actually contains two Schrodinger 
equations and the solutions of both must yield the same 
eigenvalue E 2 . 

The single kink (and antikink) backgrounds are odd 
functions of z, we see that under z — > — z, their first 
order equations transform into 



(d x ±g<j>)0 ± =TEfe- 



(13) 



That is, the parity reserved positive energy solutions are 
the parity un-reversed negative energy solutions. In other 
words, since kink and antikink are parity reversed func- 
tions of each other, the positive energy solutions on the 
kink are the negative energy solutions on the antikink; 
the negative energy solutions on the kink are the pos- 
itive energy solutions on the antikink. Further, since 
the derivative of an odd function is an even function 
we observe that the corresponding Schrodinger equation, 
Eq. (|11[) . is invariant under parity transformation: hence, 
if the energy eigenstatcs turn out to be non-degenerate 
(they are, as wc will sec below), they must be of a definite 
parity. 



For even cf>, the first order equations ([7]), (JU) transform 
under parity z — ► — z into 



{d z T9<l>)0±=±E0 T . 



(14) 



and hence 0+(z) = f3-(—z). This includes the case of the 
kink-antikink background. An alternate way to see this 
is that d z (j) is an odd function of z, and the Schrodinger 
equation for (3-(z) is identical to that for /?+(— z). Hence 
if we have a solution to Eq. (TTTj) for /3+(z) for the kink- 
antikink background, (3-{z) = f3 + {—z) will be a solution 
for the /?_ Schrodinger equation with the same value of 
E 2 . In what follows, for the kink-antikink background, 
we will simply work with the /3+ equation. 



III. FERMION BOUND STATES ON A KINK 

We begin by solving the Schrodinger equation for a 
fermion on a single kink. 



where 



d 2 J± + V K , ± {z)l3 ± =E 2 f3± 



V K ,±(z) ee g 2 -g(g± l)sech 2 



(15) 



(16) 



For any value of g > 0, Vk,+ has the shape of a poten- 
tial well with asymptotic maximum of g 2 , and minimum 
value of —g at z = 0. We know from quantum mechan- 
ics in 1 dimension that every non-positive potential that 
tends to zero asymptotically necessarily has at least one 
bound state. Hence Vk,+{z) has at least one bound state 
for every g. Also, since Vk,+ (z) gets deeper with increas- 
ing g, we expect more and more bound states to appear 
with larger values of g. This expectation will be con- 
firmed below. However, we also need a non-trivial bound 
state of the (3- Schrodinger equation which has the same 
energy eigenvalue as for (3 + . Only then will /3± solve the 
first order equations, Eq. ((5J), except if E = for then we 
can take j3- = 0. For < g < 1, Vk,- is in the shape of 
a potential barrier and clearly has no bound states. This 
shows that for < g < 1, the only possible bound state 
is with E = and /?_ = 0; the solution is 



sech 9 ^ 



(17) 



More bound states do appear for g > 1 as we now find 
by explicit calculation. 

Employing the prescription in Refs. 0, 0] we write 



0± =Af±sech b zF±(z) 



(18) 



with b 2 = g 2 — E 2 7 or b = + \J g 2 — E 2 , the positive 
choice of sign to ensure square integrability. Next we 
switch variables to 



1 



(1 — tanh z) 



(19) 
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and obtain the hypergeometric equation, 

u(u-l)F£(u) + {b+l){2u-l)F^(u) + 

(b(b+l)-g(g±l))F ± (u) =0 (20) 

It can be inferred that the arguments of the hypergeo- 
metric function F[a±, /3±; 7±; u] must be 



a± 



= b 



p± = b 
7± = b 



1 ( 1 

1 

2 
1 



1 

g± - 
y 2 



(21) 



Observe that the (g ± 1/2) actually comes from tak- 
ing a square root, so it ought to be contained within 
an absolute value sign, \g ± 1/2|; but including a and 
/3 without the absolute value sign already covers both 
cases g ±1/2 > and g ±1/2 < 0, since the hyperge- 
ometric function obeys the symmetry F [a±, /3±; 7; u] = 

The general solutions for (3± are therefore 



As z 



0±( z ) = Cisech b z F[a±, 0±; j±; u] 
-C 2 e bz F[a ± - 7 + 1, 0± - 7 + 1; 2 - 7; u] 

-> +00, tanhz - 



(22) 



-1 and from Eq. (jAlj) the hy 



pergeometric function after the e term goes to 1. As 
a result, we see that the second C2 term becomes un- 
bounded because of the e bz factor. Hence we need to set 
C2 = for normalizability. 

As z — > —00. we use the identity in Eq. (|A2[) to inform 
us that, 



lim (3 + {z) = Af+ [e 



,62 



r[6+i]r[-6] 



-bz 



r[s- 
r[6 



m-g] 

l]T[b] 



r[b + g + l]T[b-g] 



(23) 



lim p_(z) = Af- [e 



-bz 



. r[b±m-b] 
r[s]r[i - g] 
r[b + i]r[6] 

r[b±g}T[b-g±l] 



(24) 



The e~ bz term would be unbounded if its coefficient is 
finite. Recalling that the gamma function has poles at 
the negative integers and zero, we can then set the e~ bz 
term to zero by requiring that the argument of one of 
the gamma functions in the denominator to be a negative 
integer or zero. Since both b + g and b + g+1 are strictly 
positive, we need 



l± 1 1 

K - g + o t r = -n± e 



(25) 



which implies 



E n+ = \Jn+(2g - n+) 
E n _ = V(n_ + l)(2 5 -(n_ + l)) 



The solution for /3± is valid only if their energy eigenval- 
ues coincide, we get the additional requirement 



n + — n- = +1 



(26) 



The range of n + is determined by noting that 6+ = g — 
n+ from Eq. (|25| and normalizability requires 6+ > 0. 
Therefore 



< n + < g 



(27) 



We then need to determine the relationship between 
the normalization constants Af± of these /3+ and /3_ so- 
lutions by plugging them back into our first order equa- 
tions © . With some algebra involving the hypergeomet- 
ric function identities (|A3[) and (|A4[) . we can verify that 
our solutions do satisfy the first order equation provided 
we have 



Afi n) 



A/ 



(») 



n 



(28) 



where n = n + labels the n th mode. 

To summarize, on the kink background the positive 
energy fermionic bound states are given by 



l3 { +\z) = -N n E n sech 9 - n z 

1 



F 



-n,2g 



1; - (1 - tanh(z)) 



0( z ) = J\f n n sech 9 ~ n z 



F 



-n + 1, 2g — n; g — n + 1; — (1 — tanh(z)) 



(29) 



(30) 



E„ = y/n(2g-n), < n < 3, n G Z + 

where wc highlight that, because — n and — n + 1 are 
negative integers or zero, we sec from (|A1I) the hyperge- 
ometric functions are really finite order polynomials in 
11= (1 — tanh z)/2. 

F [—71, 2g — n + 1; g — n + 1; u] = 

(-n) m (2g+ 1 -n) m ^ m 



m=0 



77l!((7 — 71 + 1), 



F [— n + 1, 2g — n; g — n + 1; u] = 

E 1 (— n + l) m {2g - n) r 
r. 7^ 1 



m=0 



m\(g — n + 1), 



As an example, we can recover the bound state found 
in Ref. @ by setting n = 1. 

P ( +\z) = -Af v / 2g- 1 scch ff_1 z tanhz 
0(z) = Afscch^z 

E x = ^2g~l (31) 
where M is a normalization factor. 
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FIG. 1: Kink-antikink potentials V K j? for g — 0.5 and g — 1.3. 



IV. BOUND STATES ON KINK-ANTIKINK 

As discussed below Eq. (fT4"|) . at the end of Sec. HH it 
is sufficient to find the solution for (3+(z) in the kink- 
antikink background and then set (3-(z) = (3+(—z). So 
we will only focus on finding (3 + . 

On inserting the kink-antikink background of Eq. ([3]), 
the Schrodingcr equation (fTTj) becomes 



h kW /3+ = (-a 2 + y KK ) P+ = Elp + 



(32) 



where the potentials are 



^kk = v *.+ + V K> - - g 2 

+ 2. 9 2 e^ 2L sech(z + L) sech(z - L) (33) 

where the expressions for Vk.± ar e given in Eq. (jTHJ) . The 
shape of this potential is illustrated in Fig. [1] for g = 0.5 
and 1.3. 



A. Proof of existence of bound states 

There is a theorem by Simon [8j which states that a 
potential eV(z) admits at least one bound state for all 
e > if and only if J^° oo V(z)dz < 0. 1 Applying this 
criterion to our potentials (shifted by — g 2 ), 



-t-oo 



VI KK) (z) -g 2 )dz = -4g z + 8g 2 L— 



-,-2L 



sinh(2L) 

(34) 

At large L, 8g 2 Le 2L / sinh(2L) is small compared to 4g 2 , 
and hence the integral is negative. Solving for the zero 
of the right hand side amounts to solving 



AL + 1 = e 



-IL 



(35) 



1 An elementary proof by computing the expectation value of the 
Hamiltonian with respect to some trial wavefunction, can be 
found in 0], 



But y = 4i + 1 is the tangent line to y = e at L = 0. 
That is, the only solution to the above equation, and 
hence the only instance the integral of the potential be- 
comes non-negative, is when L = 0. For all L > 0, there- 
fore, we see that the kink-antikink background, as speci- 
fied by Eq. Q . supports at least one fermion bound state 
for all non-zero values of the coupling g. Contrary to the 
claim by Postma and Hartmann [6| , we see that spin does 
not pose any obstacle to the existence of fermion bound 
states on the kink-antikink. 



B. A lowest energy upper bound 

As mentioned in the fermion zero mode (E = 0) 
solution on the kink-antikink is not normalizable, as can 
be verified by integrating directly. That means Eq is 
strictly positive. From the variational principle in quan- 
tum mechanics, we also know that the ground state en- 
ergy Eq is always less than or equal to the expectation 
value of the Hamiltonian H K -^ with respect to an arbi- 
trary square intcgrable wavefunction \ip), namely. 



E 2 < (i>\H K K\i>) 



Motivated by the fact that 

(p(z) = sech s (z + L) 



(36) 



(37) 



is the /3+ zero mode solution to a single kink at z = — L 
and the only normalizable (3+ solution to the antikink at 
z = +L is zero, we shall use tp as our trial wavefunction. 

Inserting the Hamiltonian in Eq. (|36[) and using the 
equation obeyed by the zero mode state (Eq. (fl~5|) with 
E = 0) we get 



< El < 



dz sech 2s z + sech 



— g(g — l)sech z_ + 2g 2 e _2i sech z + 



(38) 



where we have denoted z± = z ± L and also used the 
result M, [Til 



sech 2ff zdz 



r\9 + \] 



(39) 



The second term in the bracket in Eq. (|38|) gives a con- 
tribution proportional to 

2g 2 e- 2L J dz sech 29+1 z + sech z_ 

J dz e z sech 29+1 z (40) 



< 8q 2 e~ 4L 



where we have used the inequality sech z_ < 2e z ~ . The 
first term in the bracket also gives a contribution propor- 
tional to e~ iL for g > 1. However, for < g < 1, the 
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contribution is estimated using 

g(l — g) J dz scch 2ff z + scch 2 z_ 

< g)2 2g e- igL J dze 2gz scch 2 z (41) 

where we have used the inequality sech 2£, z+ < 2 2g e 29Z+ . 
The end result is 

o < e 2 < e -4L r b + i] 8g2 r dz e - sech 2 9 +i z (42) 



ln|E»-Vn(2g-n)| 



if g > 1, and 

< E 2 < e- 4gL r ^ + r ^ 9 (l - 9 )2 2 ^ / e 2 ^scch 2 z 

(43) 

if < g < 1 in the large L limit where the first term in 
Eq. I|38p dominates over the second term. 

These results provide an upper bound for the energy 
of the ground state in the kink-antikink background, the 
existence of which we proved in the previous subsection. 
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FIG. 2: Ground state and e xcited ener gy levels of fermion on 
kink-antikink near E„ = \Jn(2g — n), which are the energy 
levels on the single kink, for g = ir. Here we plot the absolute 
value of the deviation from E n to show, for the first three 
levels, the roughly exponential dependence on L, i.e. 8E n = 
\E n — E n \ ~ e~ aL , with a > 0. From dark to light, the dots 
are for n = 0, 1 and 2, with best-fit slopes of -6.28, -3.41 and 
-2.25 respectively. 



C. Numerical Solutions 

We proceed to numerically solve the fermion bound 
state on the kink-antikink. 

First we note that it is impossible for (3± to both vanish 
at the same z. Recall that first order equations arc solved 
uniquely by specifying one boundary condition for each 
[3. So if it were the case that (3+(zq) = (3-(zo) = for 
some zq, then looking at ([5]), the unique solution is simply 
/3+(z) = fi-(z) — Vz. In particular, we cannot have 
both @± go to zero at z = 0. As discussed earlier, since 
fl+{z) = j9— (— z) for the kink-antikink we can thus set 
fl±( z — 0) = 1 and rescale the solutions later if necessary. 

The eigenvalues are written as Eq = y/\2g — 1\S and, 
for 77. > 1, E n = E n (l + 5), with E n = y / n(2g - n). 
They are searched for by solving (JSJ) repeatedly with var- 
ious values of S, and watching the large \z\ asymptotic 
behavior of the solutions, as in the "shooting method" . 
All of them eventually blow up, but as one tunes S, the 
/?+ may say switch from going to negative infinity to go- 
ing to positive infinity, as z — > — oo. The exact eigenvalue 
lies between these two values of <5 where this transition 
takes place, and the search for the eigenvalue primarily 
involves narrowing the gap between these two 8s until the 
desired accuracy is achieved. 

We selected g = tt and investigated h ow the en- 
ergy levels near those of the single kink, \J n(2n — n), 
n G {0, 1, 2}, are varied as the kink-antikink separation is 
altered from L = 2.5 thru L = 7. Referring to Fig. [21 one 
can infer that the first three energy levels roughly have 
an exponential dependence on the kink-antikink distance: 
E n ~ e~ aL , for some a > dependent on n. This indi- 
cates the {E n } approach that of their single kink coun- 
terparts as L is increased, in accordance with physical 
intuition. 
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FIG. 3: Ground state energy vs. g, the Yukawa coupling, for 
L = 5. We see that E ~ e" 8 ' 369 . 



For L = 5, we varied the coupling g from 0.1 thru 4 
to examine the effect on the ground state energy eigen- 
values. Fig. [3] provides evidence that the energies de- 
crease roughly exponentially with increasing strength of 
the coupling. 

The remaining figure, Fig. [4] shows the numerical /3+ 
solution to the kink-antikink system for the ground state 
of {g, L} = {0.1, 5}. It is compared against the cor- 
responding analytic solution (3+(z) = sech 9 (z + L) for 
the single kink at z = —L; the /3+ solution for the sin- 
gle antikink at z = +L is zero. The numerical solution 
is normalized so that its approximate peak at z = — L 
coincides with that of sech 9 (z + L). 
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FIG. 4: Ground state of fermion on kink-antikink with g = 
0.1, L = 5, and £?o ~ 0.04. The solid line is the ground 
state (3+ = sech s (z + L) solution on a single kink centered 
at z — —L. The dashed line is the numerical solution to the 
kink-antikink system. 



V. CONCLUSIONS 

We have tackled the problem of solving for bound 
states of the Dirac equation in (1+1) dimensions on kink 
and kink-antikink backgrounds. The resulting coupled 
first order equations can in turn be uncoupled to yield 
two Schrodinger equations, which we solve exactly for 
the single kink and antikink case. We find that the num- 
ber of positive energy bound states on a kink is given by 
the smallest integer less than g = 2m//ra s . For fcrmions 
on a kink-antikink, we used the Schrodinger equations 
and results from non-rclativistic quantum mechanics to 
prove that at least one bound state has to exist, for all 
non-zero values of the Yukawa coupling g. We then de- 
rived an upper bound for the lowest energy squared Eq 
value which allowed us to prove that the ground state en- 
ergy of the fermion on the kink-antikink tends to zero as 
the kink-antikink separation tends to infinity (L — > oo). 
Appropriate boundary conditions for the first order equa- 
tions were devised and employed to solve numerically the 
energy eigenvalues and eigenfunctions. For the specific 
examples we looked at, the lower lying bound states ap- 
proached that of their single kink counterparts exponen- 
tially quickly as the kink-antikink distance was increased. 
Similarly, the ground state energy approached zero ex- 
ponentially quickly as one increased the strength of the 
Yukawa coupling. 

We expect our results to be valid also for the case of 
vortex-antivortex pairs, and for the case of loops of cos- 
mic string. The lowest non-negative energy state on a 
loop of cosmic string will have positive energy that is 
suppressed by exp(— cR/w) where R is the radius of the 
loop and w is a width associated with the string and c 
is a numerical constant of order unity. In cosmological 
applications, this is an enormous suppression and we ex- 
pect the picture derived on the assumption of exact zero 
modes to still hold true. Exceptions could occur if a 



loop shrinks and becomes small, or where a cusp occurs 
on a loop. For the case of superconducting strings 
the small but non-zero energy of the lowest positive en- 
ergy state means that charge carriers now have to jump 
from the Dirac sea to positive energy, requiring 2m en- 
ergy, where m is the mass of the lowest positive energy 
state. An applied electric field with strength < m 2 /e 
along the string can cause this jump as in Schwingcr pair 
production but the process is due to tunneling and is ex- 
ponentially suppressed [TB|. At stronger electric fields, 
the process would be unsuppressed. The critical value 
of the electric field for unsuppressed pair production is 
~ m 2 f exp(—cL/w)/e where e is the electric charge of the 
fermion. 

Another setting where fermion zero modes are believed 
to play an important role is in branc cosmology where 
fermions are trapped on 3+1 dimensional brancs in a 
higher dimensional bulk universe. If the fermions have 
zero modes in the brane background, it corresponds to 
massless fermions that are trapped on the brane and this 
is a possible explanation for massless standard model 
fcrmions living in a 3 dimensional space. In light of our 
results, if the brane can be thought of as a domain wall, in 
addition to the fermion zero modes, we may also expect 
other bound states to exist for a range of parameters. 
If the brane is closed or the bulk contains neighboring 
antibranes, the fermion zero modes will become bound 
states with an exponentially small mass. This may ei- 
ther be viewed as an undesirable feature of the particu- 
lar brane system, or else may be viewed as a means to 
probe brane configurations in the bulk via the properties 
of standard model fermions. 
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APPENDIX A: HYPERGEOMETRIC FUNCTION 
IDENTITIES 



In this appendix we collect various hypergeometric 
identities 0, O EH, EH used in this paper. 



F[a, 0; 7; u] = ^ 



m!(7), 



(a) m = (a){a + 1) . . . (a + m - 1), |u| < 1 (Al) 



7 



F[a,0;T,u] = 

r[ 7 ]r[ 7 - a -/3] 

r[y-a]T[y-0\ 
+(l-u) 



F[a, /?; 1 + a + /3 - 7; 1 - it] 



7 -a- g r[ 7 ]r[a + /3-7] 



r[a]r[/3] 

xF[-f - a, 7 - /3; 1 - a - (3 + 7; 1 - it], 
|arg[u]| < 7r, |arg[l - u]| < tt, 
a + /3- 7 ^0,±1,±2,... (A2) 



u ^ _i? [ a ' ^ ; T, u \ = a ( F l a + !> Z 3 ; 7; u] ~ ^[", 7; «]) 

(A3) 



(a + l-/3)(l-u)F[a + l,P; T ,u] = 
(a + 1 - 7)F[a,/3;7;«] 

+( 7 -^[a+l,/3-l; 7 ;w] (A4) 
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